Anistropic Stars: Exact Solutions and Stability 
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I report on recent work concerning the existence and stability of self-gravitating spheres with 
anisotropic pressure. After presenting new exact solutions, Chandrasekhar's variational formalism 
for radial perturbations is generalized to anisotropic objects and applied to investigate their stability. 
It is shown that anisotropy can not only support stars of mass M and radius R with 2M/R > 8/9 
and arbitrarily large surface redshifts, but that stable configurations exist for values of the adiabatic 
index 7 smaller than the corresponding isotropic value. 



I. INTRODUCTION 



A common assumption in the study of stellar structure and evolution is that the interior of a star can be modeled 
as a perfect fluid Q. This perfect fluid model necessarily requires the pressure in the interior of a star to be isotropic. 
This approach has been used extensively in the study of polytropes, including white dwarfs, and of compact objects 
such as neutron stars 0. However, theoretical advances in the last decades indicate that deviations from local isotropy 
in the pressure, in particular at very high densities, may play an important role in determining stellar properties |3|]. 

The physical situations where anisotropic pressure may be relevant are very diverse. By anisotropic pressure we 
mean that the radial component of the pressure, p r (r), differs from the angular components, pe{r) = Pip(r) = Pt(r). 
(That pe{r) — p v {t) is a direct consequence of spherical symmetry.) Of course, spherical symmetry demands both to 
be strictly a function of the radial coordinate. Boson stars, hypothetical self-gravitating compact objects resulting 
from the coupling of a complex scalar field to gravity, are systems where anisotropic pressure occurs naturally Q. 
In the interior of neutron stars pions may condense. It has been shown that due to the geometry of the tt~ modes, 
anisotropic distributions of pressure could be considered to describe a pion condensed phase configuration Q . The 
existence of solid cores and type P superfluidity may also lead to departures from isotropy within the neutron star 
interior 0]. 

Since we still do not have a detailed microscopic formulation of the possible anisotropic stresses emerging in these and 
other contexts, we take the general approach of finding several exact solutions representing different physical situations, 
modeled by ansatze for the anisotropy factor, p t — p r . Previous studies have found some exact solutions, assuming 
certain relations for the anisotropy factor [(J. Our goal here is two fold: first, to find new exact solutions which may 
better model realistic situations and explore their physical properties Q ; second, to investigate their stability against 
small radial perturbations For this, we generalize Chandrasekhar's well-known variational approach to anisotropic 
objects. We find that not only interesting exact solutions can be found 0, but that they may have a wider stability 
range when contrasted with their isotropic counterparts . 



II. RELATIVISTIC SELF-GRAVITATING SPHERES: BASICS 



We consider a static equilibrium distribution of matter which is spherically symmetric. In Schwarzschild coordinates 
the metric can be written as 



ds 2 = e v dt 2 - e x dr 2 - r 2 d6 2 - 



r sin 



(1) 
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where all functions depend only on the radial coordinate r. The most general energy-momentum tensor compatible 
with spherical symmetry is 

Ty = diagO, -p r , -p t , -p t ) . (2) 

We see that isotropy is not required by spherical symmetry; it is an added condition. The Einstein field equations for 
this spacetime geometry and matter distribution are 



-A 



"T + 4f ] - 3 = 87rp r ; (3) 



A' 1 \ 1 



with 



eX [ - - -i ) + -i = 8n P • ( 5 ) 

Note that this is a system of 3 equations with 5 unknowns. Consequently, it is necessary to specify two equations of 
state, such as p r = p r {p) and pt — Pt{p)- 

It is often useful to transform the above equations into a form where the hydrodynamical properties of the system 
are more evident. For systems with isotropic pressure, this formulation results in the Tolman-Oppenheimer-Volkov 
(TOV) equation. The generalized TOV equation, including anisotropy, is 

-T 1 = -{P + Pr)ir + ~(pt -Pr) , (6) 
ar 2 r 

1 , = m(r) + 47rrV 

2 V r{r-2m) ' 1 ' 

and 

m(r) = / 4nr 2 pdr . (8) 

Taking r = R in the above expression gives us the Schwarzschild mass, M. [This implicitly assumes that p = for 
r> R] 

In order to solve the above equations we must impose appropriate boundary conditions. We require that the solution 
be regular at the origin. This imposes the condition that m(r) — > as r — > 0. If p r is finite at the origin then v' — > 
as r — > 0. The gradient dp r /dr will be finite at r = if (p t ~ p r ) vanishes at least as rapidly as r when r — > 0. This 
will be the case in all scenarios examined here. 

The radius of the star is determined by the condition p r (R) — 0. It is not necessary for pt(R) to vanish at the 
surface. But it is reasonable to assume that all physically interesting solutions will have p r ,Pt > for r < R. 



III. EXACT SOLUTIONS 

In ref. , we obtained several exact solutions for different forms of the pressure anisotropy. Our solutions fall into 
two classes: i) p = constant, and ii) p oc 1/r 2 . Given the limited space, we will restrict the presentation to the latter 
case. Interested readers should consult ref. for details. 

Consider stars with energy density modeled as 

where both a and b are constant. The choice of the values for a and b is dictated by the physical configuration under 
consideration. For example, a = 3/7 and 6 = 0, corresponds to a relativistic Fermi gas, as in the Misner-Zapolsky 
solution for ultradense cores of neutron stars 0- If we take a = 3/7 and b =/= then we have a relativistic Fermi gas 
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core immersed in a constant density background. For large r the constant density term dominates (r 2 3> a /3b), and 
can be thought of as modeling a shell surrounding the core. We also take the pressure anisotropy to be 



Pt ~ Pr = 



(10) 



with c and d constant. 

We found it convenient to seek solutions for the metric function v(r) directly, rather than solving the generalized 
TOV equation. We then use the known functions A(r) and i/(r) to find the radial and tangential pressures. From eqs. 
012 and we find 



{u'f 



e — v | — + — 

4 2r 



A' 



1 



2 4 

Introducing a new variable y = e^, eq. II II becomes. 



e- x - 



1 *L 

.2 + 2r 



A' 



1 



(y"K A - y' { 2- + - ) e 

Since e~ A = 1 — 2m(r)/r, using eq. [j|]we find 



—A 



\2 + 2r ' e r 2 



= 8ir(p t -p r 



8iry(p t -p r )- 



= l-a-br 2 =I 2 (r), 



(11) 



(12) 



(13) 



where we defined the function Ij*(x) = 1 — a — bar. When 6 = 0, we write Jq = 1 — a. Using the expression for e in 
eq. IT2land substituting for the pressure anisotropy we find 



[br* - (1 - a)r 2 } y" + (1 - a)ry' - (a - c - dr 2 )y = 



(14) 



The full solution of eq. ^] with a, b, c, d ^ is in 0- Here, we will only consider solutions with b = d = 0. In this 
case, the total mass is M = aR/2, and exp[— A] = 1 — a. Since for any static spherically-symmetric configuration we 
expect (2M/ R) cr it < 1, we must have a < 1. (Also, the metric coefficient g rr becomes infinite when a = 1). 

Since we want to construct stars with finite radii and density in the context of anisotropic pressure, we impose 
boundary conditions such that p r (R) = 0. With b = d = 0, eq. ^] reduces to an Euler-Cauchy equation, 



(1 - a)r 2 y" - (1 - a)ry' + (a - c)y = 0. 



(15) 



The solutions of this equation fall into three classes, depending on the value of q (q is real, q = 0, or q is imaginary), 
where 



_ (l + c-2a)5 

We will only show results for q real. (See ref. for the other cases.) 
The solution for y is 



(16) 



/ r \ 1 +9 

9 = A+ U 



(17) 



with the constants A + and A_ fixed by boundary conditions. For the case under consideration here (b = d = 0), the 
boundary conditions are 



dr R 



Applying the boundary conditions we find 



l-3/ 2 
4ql 



and A- 



A+(q^-q). 



The radial pressure for this case, after substituting the expressions for A+ and A~, is 



8irp r 



(3I 2 -l) 2 -4q 2 lS 



R 2 i 



■2,i 



(3I 2 - 1 + 2qI 2 )R 2 i + (1 - 3I 2 + 2ql 2 )r 2 i 



(18) 



(19) 



(20) 
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FIG. 1: Radial pressure as a function of r for p and (pt — p r ) oc r and q real. 

We note that the boundary conditions automatically guarantee that p r {R) = 0. The radial pressure is always greater 
than zero provided a < 2/3 and a 2 > 4c(l — a). Since by definition a > 0, the second condition implies c > 0. Thus, 
this model does not allow for negative anisotropy. Further, since we are considering the case q > 0, we must impose 
the condition 1 + c < 2a. Combining the two inequalities for a and c, we obtain, 2a — 1 < c < a 2 /4(l — a). Since 
we have < a < 2/3 we find that < c < 1/3. We note that for the anisotropic case the maximum value of a is 
2/3, corresponding to a 33% increase when compared with the isotropic case (a = 3/7). In figure 1 we plot the radial 
pressure, p r , as a function of the radial coordinate r, for a — 3/7 and several values of c. Note that for this choice of 
a, the inequality c < a 2 /4(l — a) imposes that c < 0.08 for positive pressure solutions. This can be seen in the figure. 
For larger anisotropies, no static self-gravitating stable configuration is possible. For r <^ R, we find 

Sttt-V = 31$ -1- 2qll . (21) 

Choosing a — 3/7 we recover, in the limit c — > 0, the Misner-Zapolsky solution 0, with p r = l/(567rr 2 ) = p/'i. 

IV. STABILITY 

For sake of brevity, we will skip most details of our generalization of Chandrasekhar's formalism to 
anisotropic spheres. Readers can consult ref. [U for details. As in Chandrasekhar's original for- 
malism, we limit our study to small radial and baryon- number conserving perturbations. Writing 
p = p + Sp, p r — p ro + Sp r , p t — pt Q + 5pt, A = A + <5A, and v = v a + 5v we find that, 
to first order in v = d^/dt (£ is the Lagrangian displacement) and using the zeroth-order equations, the perturbation 
in the radial pressure satisfies, 

S P r = -Prj - IVrf^- (VW 2 £)' + yH ^ , (22) 

with 7 being the adiabatic exponent defined as 

1 

^ Pr(dn/dp r ) 

and II = pt~p r - We now assume that all perturbations have a time dependence of the form e lut . Further, considering 
5\, 5v, 5p, 5p r and SH to now represent the amplitude of the various perturbations with the same time dependence 
we obtain, after using the zeroth-order equations, the pulsation equation governing the radial stability of anisotropic 



(f>- 



dn 
Pr)^— 

Op r 



(23) 
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a = 2/9 


lu 2 R 2 = 0.95(7 -1-79) + (101.1 -52.6 7 )c 


a = 2/7 


lu 2 R 2 = 2.3(7 - 1-83) + (122.3 - 59.3 7)0 


a = 3/7 


lu 2 R 2 = 0.57(7 - 1-93) + (15.2 - 5.17)0 


a = 3.4/7 


lu 2 R 2 = 0.4(7 - 2.6) +(8.9 - 2.3 7 )c 


a = 3.49/7 


lu 2 R 2 = 0.36(7 - 2.76) + (8.0 -1.97 7)0 



TABLE I: u> 2 vs. c for given values of a. 



a = 2/9 


C,nax = 0.0016 


7 C = 1.79 -6.87 c 


a = 2/7 


c max = 0.0028 


7 C = 1.83 - 13.39 c 


a = 3/7 


Cmax = 0.083 


7c = 1.93 - 5.55 c 


a = 3.4/7 


Cmax 0.11 


7c = 2.6 - 2.84 c 


a = 3.47/7 


Cmax — 0.12 


7c = 2.75 - 7.29 c 



TABLE II: 7c vs c for given values of a. 



stars 



r 

1 



e (A +3, )/2 7 ^ (r 2 e -, /2 e y 



(24) 



8Tre*°(n o +p ro ){p +p ro )£ 



(Po+Pro) 



KP°+Pr ) r 2 (p + Pro ) 



-{\ +2v )/2 



,(Ao+2, )/22 eno / ^ 
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The boundary conditions imposed on this equation are 



£ = at r = and (5p r = at r = R. 



(25) 



The pulsation eq. I|24() . together with the boundary conditions eq. (|25|) . reduce to an eigenvalue problem for the 
frequency u and amplitude £. One multiplies this equation by r 2 £e^ A+l/ - ) / 2 and integrates over the entire range of r, 
using the orthogonality condition 



(26) 



where £ l and £ J are the proper solutions belonging to different eigenvalues w 2 . 

We now apply this equation to the exact solutions with p oc 1/r 2 described above. Using the trial function 
£ = r 2 (p + p r }e v we found that all integrals could be computed exactly. In table 1 we present results for the 
frequencies of radial oscillations u> 2 as a function of the anisotropy parameter, c, for given values of the density 
parameter. Instabilities set in for uj 2 < 0. This can be used to find the critical value for the adiabatic index, 7 c (c), 
and the maximum value for the anisotropy paramenter, c max . We also give, in table 2, the values of "f c above which 
stable oscillations are possible. Here we see that the effect of a positive anisotropy is to reduce the value of 7, thus 
giving rise to a more stable configuration when compared with the corresponding isotropic model. In particular, for 
the Misner-Zapolsky solution (a = 3/7), we find that a small positive pressure anisotropy in the equation of state 
improves the neutron star's core stability. 



V. CONCLUSIONS 



I have presented a summary of results obtained with Krsna Dev concerning the existence of self-gravitating spheres 
in General Relativity with anisotropic equations of state, aka anisotropic stars. I have also presented a summary of 
our investigation of their stability, based on the extension of Chandrasekhar's celebrated variational formalism for 
isotropic spheres to those with anisotropic energy-momentum tensors. 
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Although the search for exact solutions restricts the forms of anisotropy we could treat, our results illustrate the fact 
that, indeed, pressure anisotropy may greatly affect the physical structure of the star, leading to several observational 
effects. Most importantly, the absolute stability bound 2M/R < 8/9 can be violated, and the star's surface redshift 
may be arbitrarily large (z s > 2). It is thus conceivable that objects which are observed at large redshift may actually 
be closer than we think due to anisotropic distortions. They may also be more stable then we think, especially if 
pressure anisotropy exists near the stellar core. 

Here, I have presented results only for one of the cases we treated, stars with energy density scaling as 1/r 2 . These 
are of interest as they include ultra-relativistic equations of state used to model the cores of neutron stars. Perhaps 
the most important lesson of our study is that one must keep an open mind as to whether isotropy is or not a justified 
assumption to describe stellar matter. Until we have a better microscopic understanding of what truly goes on inside 
ultra-dense compact objects, isotropy should be taken with a grain of salt. Especially if some of these objects contain 
bosonic condensates at their cores, as is the case for several models of neutron star interiors and for a whole class of 
fully-anisotropic hypothetical objects known as boson starsQ. 
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